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ABSTRACT: 

Percolation with edge-passage probability p and first-passage percolation are studied 
for the n-cube B„ = {0, 1}" with nearest neighbor edges. For oriented and unoriented 
percolation, p — e/n and p — 1/n are the respective critical probabilities. For oriented 
first-passage percolation with i.i.d. edge-passage times having a density of 1 near the 
origin, the percolation time (time to reach the opposite corner of the cube) converges 
in probability to 1 as n ^ oo. This resolves a conjecture of David Aldous. When 
the edge-passage distribution is standard exponential, the (smaller) percolation time for 
unoriented edges is at least 0.88. 

These results are applied to Richardson's model on the (unoriented) n-cube. Richard- 
son's model, otherwise known as the contact process with no recoveries, models the spread 
of infection as a Poisson process on each edge connecting an infected node to an unin- 
fected one. It is shown that the time to cover the entire n-cube is bounded between 1.41 
and 14.05 in probability as n — > oo. 
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1 Introduction and notation 



Percolation theory, broadly speaking, is the study of connectivity in a random medium. 
Grimmett (1989) gives the following example of a question addressed by percolation 
theory. "Suppose we immerse a large porous stone in water. What is the probability 
that the centre of such a stone is wetted?" The mathematical model for this is graph- 
theoretic. Small volume elements of the stone become vertices of the integer lattice 
(of dimension three in this case). Neighboring volume elements may be connected by 
a channel broad enough to allow the flow of water, or they may not be. Model this 
by letting each pair of neighboring vertices independently have a connecting edge with 
probability p, for some parameter p. The center of the stone is then wetted if and only if 
the connected component of this random subgraph containing the center extends to the 
surface. 

Oriented percolation is a variant on this, where each edge has a particular orientation, 
and the water may pass only in that direction if at all. For example, the stone may be 
subjected to water only from above, with flow of water from one volume element to a 
neighboring one occurring (due to gravity) only when the latter is lower. First-passage 
percolation is a similar model, the difference being that the passage of water through a 
channel is not simply a yes or no event, but takes an amount of time depending on the 
breadth of the channel. Thus each edge between neighboring vertices, instead of being 
randomly included or excluded, is assigned a random (i.i.d.) passage time. The question 
is not whether, but when the center of the stone first gets wet. 

Richardson's model is a stochastic process on a graph that begins with one vertex 
infected and evolves by transmission of the infection according to an i.i.d. Poisson process 
on each edge: if the Poisson process on the edge connecting v to w has a point of 
increase at time t and one oi v or w was infected before time t, then both are infected 
after time t. Questions about Richardson's model can be reduced to questions about 
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first-passage percolation with i.i.d. exponentially distributed passage times. In Pemantle 
(1988), questions about recurrence or transience of random walks in random environments 
on trees are essentially reduced to oriented first-passage percolation on trees. Lyons 
(1990) connects random walks and electrical networks with percolation on trees. These 
equivalences (along with the desire to generalize whenever possible) are among the reasons 
to study percolation on graphs such as trees and n-cubes which do not model any actual 
stones. 

Classical percolation theory, i.e., on Z^, is geometric. The main arguments — counting 
countours, piecing together sponge crossings, viewing the process from its left edge — are 
all pictorial. Infinite binary trees were introduced as a way to get graphs which were 
in some sense limits of Z'^ as c? ^ oo. Trees are in general easier to study than integer 
lattices. For example, the first-passage percolation problem for binary trees is essentially 
a large deviation calculation: if every path of length n in the tree were disjoint then 
the passage time would be the minimum of 2" different sums of n i.i.d. random vari- 
ables; such a calculation is quantitative, requiring no picture. The argument in Pemantle 

(1988) consists mostly of showing that the slight overlapping of paths is inconsequential. 
Unfortunately, trees are locally quite different from lattices (despite being called Bethe 
lattices) and are therefore not very satisfactory in the role of limiting lattices. The n-cube 
is an alternative way to capture the high-dimensional limiting behavior of integer lattices 
without altering the local connectivity properties. The solutions presented here to per- 
colation problems on n-cubes use large deviation and second moment estimates and are 
thus closer to the tree case than to the integer lattice case. In this respect, the method is 
very similar to the arguments used by Cox and Durrett (1983) for the analogous problem 
on (Z+)'^ in high dimensions. The results resolve affirmatively conjecture G7.1 of Aldous 

(1989) , despite a remark there indicating that second moment methods shouldn't work. 
(To be fair: the second moment method bounds the probability in question away from 
but a variance reduction trick is needed to get it equal to 1.) 
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Notation will be as follows. Let Bn be the Boolean algebra of rank n whose elements 
are ordered n-tuples of O's and I's. It has a bottom element = (0, ... ,0) and a top 
element i = (1, . . . , 1). It is also useful to regard elements of Bn as subsets of {1, . . . n}, 
where the set A corresponds to the sequence with a 1 in position i if and only if i & A. 
The collection of all ^4 of a given cardinality k is called the k*^ level of Bn- The n- 
dimensional cube, or n-cube, is the graph whose vertices are the elements of Bn and 
whose edges connect each set A to ^4 U {j} for each j ^ A. Sometimes we must think of 
the edges as oriented from A to AU {j}. Representing the vertices of the n-cube as the 
standard basis in makes each such edge parallel to the unit vector e(j) connecting 
to (0, . . . , 0, 1, 0, . . . , 0) with the 1 in position j. Then a path from to 1 of length n 
can be represented as a permutation (7r(l), . . . , 7r(n)) where the kth edge in the path is 
parallel to e{7r{k)) and in fact connects the set A — {7r{i) : i < k} to the set A U {7r(A;)}. 
Similarly, a path connecting a set A to a set S D A may be represented as a permutation 
ofB\A. 

Let be i.i.d. random variables as vw ranges over all oriented edges of the n-cube. 
The first-passage time from 6 to i, or percolation time, is defined as the minimum over 
paths ^ — {0,Vi,V2, ■ ■ ■ , Vn-i, 1) from to 1 of Tn{'y), where T„(7) is the sum of Xy.^^.^^ 
along the edges of 7. In the apphcation of first-passage percolation to Richardson's 
model, the common distribution of the X^^'s is exponential with mean 1. For our basic 
first-passage percolation result itself, the assumption of exponential edge-passage times 
simplifies the derivation of certain large deviation estimates but is not necessary; we 
state the theorem for more general distributions of the passage time, though probably 
still greater generahty is possible. For the integer lattice in high dimensions, Kesten 
(1984) gives sufficient conditions on the behavior of the common density of the edge- 
passage times near the origin for similar calculations to work. It is unlikely that the 
same conditions are sufficient for our problem; at any rate, our paper does not address 
this issue. 



5 



Of course, for ordinary percolation, the common distribution of the X's takes on the 
two values "open" and "closed" with respective probabilities p and I — p, for some p. 
The basic question then is to compute the probability that there is a path from to i 
consisting only of open edges. 

The rest of the paper is organized as follows. Section 2 contains a description of the 
second moment method, and its enhancement by a variance reduction technique, that 
is used to analyze oriented percolation and oriented first-passage percolation. Section 2 
also presents lemmas that count oriented paths in the n-cube. 

Section 3 discusses oriented percolation and oriented first-passage percolation. The 
probability of oriented percolation approaches a limit when n — > oo and np is constant, 
as given in the following theorem. 

Theorem 3.2: Let each edge of Bn be independently open with probability p — c/n. 
Then P(6 is connected to i by an oriented open path) converges to a limit as n ^ oo. 
The limit is if c < e and is (1 — x{c))'^ if c> e, where x{c) is the extinction probability 
for a Poisson(c) Galton-Watson process, namely, the solution in (0, 1) to x — e'^^^"^). 

For the oriented first-passage time we have 

Theorem 3.5: Let the edges ofBn he assigned i.i.d. positive random passage times with 
common density f , and assume that \f{x) — 1| < Kx for some K and all x > 0. Then 
the oriented first-passage percolation time T — T^") for Bn converges to 1 in probability 
as n ^ oo. 

Section 4 treats unoriented percolation, for which the critical probability is shown to 
be 1/n. The result here, which may already be part of the percolation folklore, is 

Theorem 4.1: Let each edge of Bn be independently open with probability p — c/n. 
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< c < oo. Then P(0 is connected to i by an (unoriented) open path)^ (1 — x(c)y, 
where x{c) is, as in Theorem 3.2, the extinction probability for a Poisson(c) Galton- 
Watson process. 

Section 5 gives an argument of Durrett (personal communication) provideing a lower 
bound for the unoriented first-passage time when the common distribution of the passage 
times is exponential by comparing the process to a branching translation process (BTP), 
which is similar to a branching random walk. 

Theorem 5.3 (Durrett): As n oo, the time Tn of first population of 1 in BTP 
converges m probability to ln(l + ^2) = 0.88. Consequently, P(r„ < ln(l + ^2) -e) ^ 0. 

Since the oriented first-passage time is an upper bound for the unoriented first-passage 
time, the unoriented first-passage time is thereby bounded as n — >^ oo between two fairly 
close constants, namely, 0.88 and 1. We remark that for oriented first-passage times, 
lower bounds are easy (first moment calcination) while upper bounds equalling these 
lower bounds are more difficult (second moment estimates). For unoriented first-passage 
times, we do not know how to bridge the gap between the lower and upper bounds. 

Finally, the sixth section discusses the cover time in Richardson's model, i.e. the 
first time that all sites are infected. As n — > oo, the cover time is bounded in prob- 
ability by a constant; this is shown by improving Theorem 3.5 so as to get an expo- 
nentially small bound on the probability of a vertex failing to be reached by a time 
c = 3 -I- 2 In (4 + 2\/3) = 7.02. This is certainly not sharp though it improves on the 
previous best upper bound for the covering time which was of order In n. A lower bound 
in probability of ln(2) + | ln(2 + y^) is also given, showing that the covering time is 
separated from the single- vertex first-passage time of 1. 
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Theorem 6.4 and Corollary 6.3: For any e > 0, P(A(2c + e) = Bn) ^ 1 as 

n oo, where c = 3 + 21n(4 + 2a/3) = 7.02. On the other hand, for any e > 0, 
P(^(iln(2 + V5) + ln2-e) = S„) ^ as n ^ oo. 

A future paper will address the discrete-time analogues of these problems, which are 
closely related to the so-called broadcasting problem discussed by Feige, et al. (1990) and 
others. 

2 Preliminaries: the second moment method and 
path counting 

2.1 The second moment method, with variance reduction 

Aldous (1989, Lemma A15.1) gives the following lemma as the cornerstone of the so- 
called second moment method. The proof is a simple application of the Cauchy-Schwarz 
inequality. 

Lemma 2.1 (Second moment method) Let N be a nonnegative real random variable 
with EN^ < oo. Then P(7V > 0) > (ENf/EN'^. □ 

In our applications, the random variable N = Nn will be indexed by the dimension 
n of the cube under consideration. When the variance of N is o{(ENy) as n tends to 
infinity, the inequality shows that P(A'" > 0) = 1 — o(l). On the other hand, if we can 
only show the variance of N to be 0((EA^)^), then the conclusion is weaker, namely, that 
P(A'" > 0) is bounded away from 0. One of the purposes of the present work is to point 
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out that the weaker conclusion may often be strengthened to the former by a simple 
variance reduction trick. Consider an auxiliary random variable Z which absorbs most 
of the variance of A^, in the sense that the conditional second moment E(iV^ | Z) is only 
(1 + 0(1)) times the conditional squared first moment (E(iV | 2'))^ uniformly over a set of 
values of Z of probability 1 — o(l). Then uniformly for those values of -Z", P(A'^ > | 
is 1 — 0(1), and hence P(A^ > 0) is also 1 — o(l). 

To apply the second moment method to oriented percolation on the n-cube, let 
be the number of paths of length n from to 1 consisting entirely of open edges. Then 
6 is connected to 1 by an oriented path of open edges if and only if > 0. For large 
enough values of p = P(edge is open), the (unenhanced) second moment method will 
get P(A^ > 0) bounded away from by showing that (EA^^)/(EA^)^ is bounded, in 
the following manner. A similar argument, to be detailed later, using the enhanced 
second moment method will yield the exact limiting value of the percolation probability 
P(A^ > 0). 

Since EA^ is the sum over paths 7 of the probability that 7 consists entirely of open 
edges and EA^^ is the sum over pairs of paths (7, 7') of the probabihty that both paths 

consist entirely of open edges, the quotient EA^^/ (EA^)^ would be precisely 1 if the events 
that 7 consists entirely of open edges and that 7' consists entirely of open edges were 
always independent. Of course they are not independent when 7 = 7', but also they are 
not independent when 7 and 7' have any edges in common, and the covariance of their 
indicators is greater the more edges that 7 and 7' share. Thus the argument rests on 
showing that the number of pairs (7, 7') sharing a lot of edges is small. 

The oriented first-passage percolation problem is handled similarly, with A^ being the 
number of paths whose total passage time is at most M; if A" > with high probabihty, 
then the passage time is less than M with high probability. 
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2.2 Counting oriented paths 



Section 2.2 contains the path-counting lemmas needed to execute the second moment 
method. By symmetry, it will only be necessary to consider the case where 7 is the path 
given by (1, . . . , n) in the permutation representation. Let f{n, k) be the number of paths 
7' that share precisely k edges with 7. Let F{n, k) = J2n>i>k fij^i 0- Most of the time it 
will suffice to bound F{n, k) and observe that /(n, k) < F{n, k). 



Lemma 2.2 Let K{n) he any function that is o{n) as n —> 00. Then f{n, k) < F{n, k) < 
(l-|-o(l))(A;-|-l)(n — A;)! as n —> 00 uniformly in k for k < K{n). Furthermore, for k > 0, 
consider paths that agree with 7 only in an initial segment and a final segment; it is these 
paths that matter for F{n, k), in the sense that all the rest of the paths only contribute 
o{{k + l)(n — /c)!), uniformly in < k < K{n). 

Lemma 2.3 Suppose k>n- n^/'^/2. Then f{n, k) < F{n, k) < {n - k + l)(2n^/^)"-^ 



Lemma 2.4 Suppose k < n — 5e(n + 3)^/^. Then, for n > 25, f{n,k) < n^{n — k)\, and, 
forn so large that \5e{n + 3f/^] < [71^/4/2], 



F{n,k) < 2n%n-k)\+\5e{n + 3f/^] (2n^/8) N"+3)^''l-i. 



The following lemma will be needed when we apply the enhanced second moment 
method. 



Lemma 2.5 Consider now a Boolean lattice Bn+2L of size n + 2L for some positive 

integer L. Let Xi and X2 be distinct vertices in level L of Bn+2L (^i^d let yi and 1/2 be 
distinct vertices in level n + L of Bn+2L- Assume that xi lies below yi for i = 1,2. Let 
H{n, L,k,xi,X2,yi,y2) be the maximum over paths 7 from xi to yi of the number of 
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paths from X2 to y2 that have at least k edges in common with 7. Let Fi{n,L,k) = 
maxj;j^a;2^j^j^2^2 i7(n, L, A;, xi, a;2, yi, 1/2)- Then for fixed L and any function K{n) that is 
o{n) as n — > OO; the following hold: 

Fi{n, L, k) = o{{k + l)(n — /c)!) uniformly in < k < K{n); (1) 

L,k)<{n-k + l)(2n^/«)('*-'=) for k > n - 71^/^2; (2) 

Fi(n, L, k) < 2n\n - k)\ + \he{n + 3)'/^] (2n'/«) for k < n - 5e(n + 3^^ 

(3) 

For the last inequality we require 5e(n + 3)^/^ < . 

The following notation is common to the proofs of Lemmas 2.2-2.4. If 7' has precisely 
k edges in common with 7 and is given by the permutation tt, then let ri, . . . , be the 
terminal positions of the shared edges. In other words, let ri be minimal so that the 
set {1, . . . , ri — 1} is equal to the set {7r(l), . . . , 7r(ri — 1)} and 7r(ri) = ri; let r2 be the 
next such value, and so on. By convention, let ro always be and r^+i always be n + 1. 
Write r = r(7') for the sequence (ro, . . . , rk+i). Write Sj = r^+i — r^, i — 0, . . . , A;. For 
any sequence Tq = (ro, . . . , r^+i) with = ro < ri < • • • < r^ < r^+i — n + 1, let C(ro) 
denote the number of paths 7' with r(7') = ro. Then it is easy to see that C(r) < G'(r), 
where 

G(r) = n(..-1)!, (4) 

1=0 

since the values 7r(rj+l), . . . , 7r(ri+Sj — 1) must be a permutation of {r^+l, . . . , rj+Sj — 1}. 
For Lemma 2.4, the more precise bound C(r) < G'i(r) will be necessary, where Gi is 
defined by 

k 

G,{r) = ll[{s,-iy.-l + 6,,s.] (5) 

i=0 

and 5i^si is 1 if Sj = 1 and otherwise. To see this inequality, recall that 7' must not 
have any common edges with 7 strictly between r^ and r^+i, and therefore, for Sj 7^ 1, at 
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least one permutation of {r^ + 1, . . . , rj+i — 1} (namely, the identity permutation) is ruled 
out for the values of n{ri + 1), . . . , 7r(rj+i — 1). Note that Sj(r(7')) can never be 2, a fact 
which is used liberally in the proofs of Lemmas 2.3 and 2.4. Of course, f{n, k) is equal 
to the sum of C(r) over all sequences r satisfying = Tq < ri < • • • < rjt < r^+i — n + 1. 
Furthermore, F{n. k) is at most the sum of G{r) over all such sequences, since G{r) 
counts the paths 7' with r C r(7'). It remains then to get bounds on these sums. We 
use the following facts about factorials. 

Proposition 2.6 (i) Factorials are log-convex, i.e., a\h\ < {a -\- j)\{h — j)\ for a>b> 

(ii) (a! - l)((a + jV- - 1) < {{a - 1)1 - l)((a + j + 1)1 - 1) as long as a > 4 or a ^ 3 
and j > 1. 

(iii) a > 6 > implies (a! - l)/{b\ - 1) > a\/b\ > (a/e)""^ 

Proof: (i) follows from (in fact, is equivalent to) the fact that (a + l)!/a! = a + 1 is 
increasing, (ii) is easy to verify. The first inequality of (iii) is trivial. To prove the second 
one, use Stirling's formula a!e~^/(^^"^ < a°'e~°■^/27^a < a\ to get 

a\/{b\{a/e)--') > ^^!^l!^^(e/a)'^-^ = (a/6)''+^/V^/(i2b) 

and taking the logarithm of the last expression gives at least 

(6 + 1/2) ln(l + 1/6) - 1/(126) > (6 + l/2)(l/6 - 1/6') - 1/(126) > 
for 6 > 2. For 6=1, the result follows directly from Stirling's formula. □ 

Proof of Lemma 2.2: For fixed r, let j = j{r) = maxj(sj — 1). Consider separately the 
cases j < n — 4k and j > n — 4k. The idea is that G{r) is small in the first case, and while 
G is not so small in the second case, there are not too many sequences with such large 
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values of j. In the first case, G{r) must be no more than (n — 4A;)!(3A;)!. To see this, note 
that by log-convexity of factorials, the product in (4) is maximized subject to j < n — Ak 
at some r for which j(r) is equal to n — 4/c, in which case since J2i{si — 1) = n — k, G{r) 
is at most [maxj(sj — — 1) — maxj(sj — 1)]! < (n — 4fc)!(3A;)!. Meanwhile, the 

number of sequences r under consideration is at most . Thus the total contribution 
from this case is at most 

{n - 4k)\{3k)\n\/{k\{n - k)\) 

n-Aky. n\ {3k)\ 



= {n-k)l 
= {n-k)l 

< {n-k)\ 



{n-k)\ {n-k)\ k\ 

[n{n - 1) • ■ ■ (n - A; + l)][3A;(3fc - 1) • • • (fc + 1)] 
(n - k){n - A; - 1) • • • (n - 4A; + 1) 

(3A;)2n 



{n-Ak + If 

The term inside the square brackets converges to as n ^ oo uniformly in A; as long 
as A; < K{n) — o{n), so for large n and k < K{n), the contribution from the case 
j(r) < n — Ak is o{{n — k)\), uniformly in these k. 

The second contribution is from the terms with j > n — 4k. Since factorials are log- 
convex, it follows (as in the preceding case) that G'(r) can be no greater than j(r)!(n — 
A; — j(r))!. The number of sequences r for which j(r) is any fixed value jo is at most 
A; -|- 1 times the number of sequences r with sq — 1 — jo, which is in turn at most 

— jo ~ l)'/((^ ~ ~ jo — So the contribution for fixed jo is at most 

{k + l)jo\in-jo-k)\{n-jo-l)\ jk + l) . „ . 

(k-mn-jo-ky. (^^°'^" - 

Changing jo to jo + 1 multiplies this by {jo + l)/{n~jQ^l) > (n — 4A'(n))/(4A'(n)), which 
tends to infinity by assumption on K{n); hence the sum is bounded by (1 + o(l)) times 
the single term with the maximum value of jo, namely, jo = n — k. The contribution 
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from this term is at most [{k + l)/{k — l)\]{n — k)\{k — 1)! = (A; + l){n — k)\. The case 
j(r) = n — k covers precisely those sequences r consisting of an initial segment and a 
final segment; thus both statements in the lemma have been proved. □ 



Proof of Lemma 2.3: Let m — n — k he the number of positions in which the edges 
of 7 and 7' differ and let j(r) this time be the number of runs of consecutive positions 
that constitute these m edges. For example, if n = 10 and k = 6 with (tq, . . . , r^+i) = 
(0,3,4,5,8,9,10,11) then j(r) = 2 because the edges not shared arc in two runs of 
consecutive positions, namely, 1,2 and 6,7. Now the inequality F{n,k) < J2rG{r) can 
be broken up as 



where g{n, k, I) is an upper bound for G{r) over sequences = ro < • • • < Vk+i — n + 1 
with j{r) = I and h{n, k, I) is the number of sequences with j(r) — I. 

As remarked before, each such run of consecutive positions contains at least two 
positions, since otherwise the single edge in the run would also be a shared edge. Thus 
I < m/2, and for fixed Z, the number of r for which j{r) — I is equal to the number of 
ways of choosing m positions out of n in Z clusters of size at least 2 each. Treating each 
cluster as a unit, there are (""7^') '^^Y^ locating the I clusters among the n — m shared 
edges. Since each cluster has at least two positions, there are m — 2l extra positions to be 
distributed among the I clusters. The number of ways to do this is the number of ways 
to drop m — 21 indistinguishable balls into I distinguishable boxes, which is (™r!r^). 



when I — [m/2 J, while the second factor is at most which is also maximized at 
I = [m/2 J. Thus, assuming for simplicity that m is even. 



F{n, k) < Y,9{n, k, l)h{n, k, I), 





The first factor is increasing in I, thus maximized 
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Using again the log-convexity of factorials, G{r) is less than or equal to the product 
of the factorials of the cluster sizes, which is maximized when all clusters have size 2 
except for one large cluster. This gives g{n, k, I) < 2^~^(m — 2/ + 2)! < 2"*/^m!, and thus 

|,(..M)M".M) < (?+i)(";7/')(„7j2'"'^™! 

< (^ + l)(8mn)'"/2. 

The same inequality X)/ fi'(^) ^! i)h{n, k, I) < (y + l)(8mn)™'/^ can be established similarly 
when m is odd. Now the hypothesis that m < implies that the bound (y + 

l)(8mn)"'/2 is at most (m + l)(2n^/8)'". □ 

Proof of Lemma 2.4: The idea this time is to define a weight w{r) with the property 
that the sum of the w{r) over all sequences = Tq < • • • < rjt+i = n + 1 is less than 1. 
Then 

^C(r) = J2w{r)[C{r)/w{r)] < max[C(r)/u)(r)], 

r r 

which will be shown to be as small as required. Let Uj denote n — Vi and let Sj = rj+i — rj 
as before, viewing the sequence of Sj's as a function of r. The weight w is defined by 

k 

i=0 



where 



(m+ l)/(m + 3) 
Q{m,s) = \ l/(m + 3) 

l/[(m + 3)(m + 2)] 



s = 1 
s = 3 
s > 4. 
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Note that Q{m, 1) + I^S^ Qijn, s) < 1 for each fixed m; hence by a simple induction 
argument the sum of w{y) over all sequences of a given length k is less than 1. It remains 
to bound for fixed n and k the quantity maXr[C(r)/w(r)]. 

The procedure will be to find for each r a quantity D{y) and a pair y^^\y^'^^ such that 
C{y)/w{y) < i(n + 3)2D(r), D{y) < Z}(rW) < {n + 3fD{Y^^^), and ^(r^^)) < D{yo). 
This will give C(r)/w(r) < |(n + 3)^D(ro), and the calculation L)(ro) = [(n — /c)! — l](n — 
/c + 3)(n — /c + 2) will follow directly from the definition of D. Then 

^ < -(n + 3)'(n - A;)! < n%n - kV. 
w{y) 2 

for n > 25, as desired. 

To begin, recall the bound (5) on C(r): 

k 

G,{y) = 1[[{s,- 1)1-1 + 6,,,^]; 

i=0 

use this together with the definition of w to get the inequality 

k 



C{y)/w{y) < G'i(r) 



YlQirii^Si) ^ 



.1=0 



Let Zi be defined to be 1 when Si — 1 and Q{ni, Si) ^ otherwise. Then Hi Q{ni, Si) ^ < 
UiiZiirii + 3)/K + 1)) < (n. Z,)[{n + 3)!/(2!(n + 1)!)] = i(n + 3)(n + 2) ^ Z^. So 

C(r)Mr)<i(n + 3fi^(r), (6) 

where D{y) :=G'i(r) nio^i- 

Now if there are i < j for which > 5 and > 4, consider the (lexicographically) first 
such pair i.j, and let 0(r) be the sequence for which Sj(0(r)) = 4, Sj(0(r)) = Sj + Sj — 4, 
and Si(0(r)) = s/ for all / 7^ i,j, where Si always refers to s^(r) unless otherwise noted. 
Then (ii) of Proposition 2.6 shows that Gi(0(r)) > Gi(r), with strict inequality if 7^ 4. 
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But also Hi Zi{4'{^)) ^ Hi Zi{r). To see this, note that each factor Zi is increasing in ni 
(even when si varies, as long as si stays at least 4 or else stays constant at 1 or constant 
at 3). Changing r to 0(r) changes the s;'s only by switching two of them that are both at 
least 4, and furthermore the values of ni for (6(r) are equal to the values of ni for r except 
when i < I < j. Since ni is increased for i < / < j, it follows that D{(f){r)) > D{r). 

Let r^-*^^ be gotten from r by iterating (j) until it is no longer possible to do so. Note 
that r*^^^ has at most one Si greater than 4 and that, by induction, D(r^^^) > D[r). 

Now if maxjSj(r(^)) = 4, then let r^^^ = r^^^; otherwise, if ^^(r^^)) > 4, let r^^^ be the 
sequence for which Sj(r(^)) = Sfe(r(-^)), Sk{r^'^^) — Si{r'^^^), and s;(r(^)) = s;(r^-'-)) for I ^ i,k. 
Clearly, G'i(r(^)) = Gi(r(-'^)). For each I ^ i,k, si is unchanged while ni is increased or 
remains the same, and so Zi is not decreased. Since Zj(r(^)) > (n + 3)~^Zj(r(^)) and 
Zfe(r(2)) > Zfe(r(^), it follows that ^(r^^)) > {n + ?>)-'^D{r^^^). Summarizing the progress 
so far, there is a sequence r^^^ with Sj(r'^^') < 4 for i < A; such that 

C{y)/w{y) < ]-{n + ^fD{v) < -(n + 3fL»(r«) < -{n + 3)^D{r^''^). 
2 2 2 

The last step is to compare D{r^'^^) to -D(ro) = [{n — k)\ — l]{n — k + 3){n — k+2), where 
ro is the sequence (0, 1, 2, . . . , A;, n + 1). Let /cs = A;3(r'^^)) be the number of i < A; for which 
Si(r(^)) = 3 and let k4 — A;4(r(^)) be the number of i < /c for which Si(r(^)) = 4. Then 
Sfc(r(2)) = n + l-k-3k4-2k3. Now the inequalities (a! - l)/(6! - 1) > a\/b\ > {a/ef-^ 
for a > 6 > from (iii) of Proposition 2.6 give 

L>(ro)/L'(r(2)) 

> [{n-k)\- l\/[h^\n + 2,f^\n + ?,f\{n -k-Zk^- 2^:3)! - 1)] 

> (5(n + 3)2/3)3fc4+2fe3 /r5fc4^^ + 3)2fc4+fc3i 
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> 1, 

where the inequahty {n — k)/e > 5(n + 3)^/^ comes from the hypothesis of the lemma. 
This proof of -D(ro) > -D(r(^)) has assumed n — k — 3k4 — 2ks ^ 1, but can easily be 
modified to treat the contrary case. 

Putting all this together gives C(r)/'u;(r) < — A;)! for all r, proving the f{n,k) 
part of the lemma. 

Finally, using Lemma 2.3, which is valid because n — 5e(n + 3)^/^j + 1 > n — 
F{n,k) < Yl fin,l) + F{n,n-\5e{n + 3y/^]+l) 

n-5e(n+3)2/3>i>fe 

< 2n^(n-k)\+\5e(n + 3)^/^] (2n^/8) -\ 

□ 

In order to prove Lemma 2.5, the following comparison is needed between overlaps of 
pairs of paths connecting to 1 and overlaps of pairs of paths connecting Xi to 

Lemma 2.7 Let xi,X2 be distinct vertices at level L of Bn+2L o,nd let yi,y2 be distinct 
vertices at level n + L of Bn+2L, with Xi below yi and X2 below 1/2 ■ Fix an oriented path 7 
connecting Xi to yi. For < k < n and i = 1,2, let Aj(7, k) be the set of oriented paths 
7' from Xi to yi that share exactly k edges with 7. For i = 1,2, let Ai [= [J^^qAi{j, k)] 
be the set of all oriented paths from Xi to yi. Then there is a bijection (f) from Ai to A2 
such that (i) the set of edges that ^(7') has in common with ^ is a subset of the edges 
that 7' has in common with 7; (ii) hence if ^' e ^1(7, A;) then 0(7') G ^2(7,^') for some 
j < k; and (iii) 2/7' and 7 share either their first or last edge, then the inclusion in (i) 
(and hence the inequality in (ii)) is strict. 
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Proof: Viewing vertices of Bn+2L as subsets of {1, . . . , n + 2L}, the path 7 is represented 
by a permutation of yi\xi. Assume without loss of generahty that xi = {1, . . . ,L} and 
yi — {1. . . . ,n + L}, and that the permutation representing 7 is in fact the increasing 
permutation (L + 1, . . . , L + n). We shall express the desired as a bijection, call it h, 
between permutations of {L + 1, . . . , L + n} and permutations of y2\x2- 

Let 7' connecting xi to yi be represented by a permutation tt = (7r(l), . . . , 7r(n)) of 
{L + 1, . . . , L + n}. Wc need to define a corresponding permutation h{7i) = {h{7T){i),i = 
1, . . . , n) of y2\x2- The idea is this: h tries to copy vr, but is required only to copy elements 
of y2\x2', so it replaces in corresponding order the elements of {L + 1, . . . , L + n} that are 
not elements of y2 \ X2 by elements of 7/2 \ ^2 that are not elements of {L + 1, . . . , L + n}. 

The construction of h can be expressed more formally as follows. Let / := {1, . . . , n}n 
[(2/2 \ X2) — L] denote the set of indices i, 1 < i < n, such that L + i e y2 \ X2; thus 
the points L + i, i E I, form the intersection of {L + 1, . . . , L + n} and y2 \ X2- There 
are m :— n — \I\ elements in {L + 1, . . . , L + n} that are not in y2 \ X2; label them in 
increasing order as /ci < A;2 < • • • < /Cm- Label the m elements that are in y2 \ X2 but not 
in {L + 1, . . . , L + n} as fc^ < < • • • < A;^. 

Let TT = . . . , 7r(n)) be a permutation of {L + 1, . . . , L + n}. If 7r(i) G y2 \ ^2, 

i.e., if 7r{i) — Lei, then let h{7r){i) — 7r{i) E L + I. If 7r{i) ^ 2/2 \ ^2, i.e., if 7r{i) — kt 
for some 1 < t < m, then let /i(7r)(i) = k'^. It is easy to see that this yields a bijection 
between permutations of {L + 1, . . . , L + n} and permutations of y2\x2- It remains to 
show that it has the required properties. 

Let 7' e Ai be represented by tt and 0(7') G A2 by h{n). Then i belongs to r(7') 
(defined as a subsequence of (1, . . . , n) in the obvious fashion) if and only if 

{7r(l), . . . , 7r(i - 1)} = {L + 1, . . . , L + i - 1} and 7r(i) = L + i. (7) 
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On the other hand, i e r(0(7')) if and only if 

{/i(7r)(l),...,/i(7r)(i-l)} = {l,...,L + i-l}\x2 and h{7r)(i) ^ L + i; (8) 

in particular, X2 C {1 . . . , L + i — 1} is necessary for (8). 

We claim that r(0(7')) C r(7'), i.e., that for each i G {l,...,n}, (8) implies (7). 
Indeed, suppose (8) holds for a given value of i. Since h{n){i) = L+i G {L + 1, . . . , L+n}, 
we have n{i) = h{Ti){i) = L + ihy our construction of hijx). Also, {/i(7r)(l), . . . , h{n){i — 
1)} is the union of some subset of {L-\-l, . . . L+i— 1} and an initial segment of {k[, . . . , k'^} 
(namely, {1, . . . , L} fl {k[, . . . , k'^}). Therefore {7r(l), . . . , 7r{i — 1)} is the union of some 
subset of {L + 1, . . . , L + i — 1} and an initial segment of {ki, . . . , km}, and so equals 
{L + 1, . . . , L + i — 1}. Thus (7) is established, and the proof of (i) is complete, (ii) 
follows immediately from (i). 

To finish the proof of the lemma, observe that since Xi and X2 are distinct, 0(7') can 
never share its first edge with 7. Similarly, yi and 7/2 are distinct, so 0(7') can never share 
its last edge with 7. Thus if 7' shares either its first or last edge with 7, the inclusion 
r(0(7')) C r(7') must be strict. □ 

Proof of Lemma 2.5: Fix Xi,yi {i — 1,2) and a path 7 connecting xi and yi. Use 
Lemma 2.7 to get a bijection from paths connecting xi and yi to paths connecting X2 
and y2 with the properties stated therein. Now the interval in Bn+2L between xi and yi 
is isomorphic to Bn- Hence the number of paths 7' connecting xi and yi and sharing 
at least k edges with 7 is just F{n,k), and Lemmas 2.2 - 2.4 may be used to bound 
this. Now since is a bijection and 0(7') shares at most as many edges with 7 as 7' 
does, this immediately gives H{n, L,k,xi,X2,yi,y2) < F{n,k)] maximizing over xi, 
X2, and 1/2 and applying Lemmas 2.3 and 2.4 estabhshes (2) and (3). To show (1), note 
that if 7" connects X2 to 1/2 and shares at least k edges with 7 then 7' = (tr^{^/') cither 
shares strictly more edges with 7 or shares exactly the same edges in which case the 
shared edges include neither the first nor last edge. The number of 7' in the former 
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category is at most A; + 1), while the number of 7' in the latter category is at most 
o{{k + l){n — k)\) uniformly in < /c < K{n) according to the last part of Lemma 2.2. But 
F{n,k + 1) < {l + o{l)){k + 2){n-k-l)\ = o{{k + l){n~ k)\) uniformly in < A; < K{n) 
according to the first part of Lemma 2.2, and the desired conclusion follows. □ 

3 Oriented percolation and oriented first-passage per- 
colation 

3.1 Oriented percolation 

The first apphcation of Lemmas 2.1 - 2.5 will be to ordinary oriented percolation. This 
means that the edges are independently open with probability p and closed otherwise. 
The problem considered here is to determine how large p should be as a function of the 
dimension n in order that and 1 may be connected with substantial probability. It 
turns out that the answer is p — e/n in the sense that ii p — c/n then for c < e the 
probability that is connected to 1 tends to as n — > 00, while for c > e, the probability 
that is connected to 1 approaches a positive limit which we calculate. (It is easy to see 
that the limit is not 1 since the disconnection probabihty is at least the chance that is 
isolated, namely, (1 — c/n)" — > e"'^.) In particular, the limiting connection probability 
is = 0.8416 at the critical value c = e. We consider oriented percolation mainly as a 
warm-up to the arguments used to analyze oriented first-passage percolation. We do not 
intend for our results to be viewed as a complete analysis of the threshold behavior of 
oriented percolation with edge probabilities c/n as c crosses the critical value e. 

The next lemma derives a lower bound of (e — l)^/e^ = 0.400 on the percolation 
probability in the critical region by using the unenhanced second moment method. The 
bound is not tight but will be a crucial ingredient to the proof of the sharper result. 
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Lemma 3.1 Fix L G M. Let each edge of Bn be independently open with probability 
p — e/{n + 2L). Then 

{e-lf 



liminf P(0 is connected to 1 by an oriented open path) > 



e2 



Proof: Observe that EiV^ is the sum over pairs of paths of the probabihty that both 
paths are open. By symmetry this is the same as n! times the sum when the first path is 
fixed, say as the path 7 corresponding to the identity permutation. The probabihty that 
both 7 and 7' are open depends only on the number k of edges they share and is equal 
to p2n-fe xhus by Lemma 2.1 it suffices to find a finite upper bound for 



(ENf ■ (p^niy 



^[/(n,fc)(^— )Vn!]. (9) 



We bound this in three pieces, corresponding to Lemmas 2.2 - 2.4; of course we may 
assume n to be as large as needed. When A; < 12 Inn, the summand in the final sum of (9) 
is by Lemma 2.2 at most (l + o(l))(A; + l)e-*^(n + 2L)'=/[n(n-l) ■■■{n-k + 1)] uniformly 
in A; as n ^ 00. This is at most (1 + o{l)){k + l)[(n + 2L)/{n — 12 lnn)]'^e~*^ and so the 
sum over A; is at most (l + o(l))(e/ [e - (n + 2L)/(n - 121nn)])^ = (1 + o(l))eV(e - 1)^ 

The contribution from any term with k > n — n^/^/2 is at most 

{n-k+ l)(2n^/S)"-'=[(n + 2L)/ef/n\ 

according to Lemma 2.3. Using the inequality n\ > (n/e)"n^/^ gives an upper bound for 
these terms of (n - A; + l)(2en-^/^^^)"-'=n-V2(l + The sum over k is then at 

most [1/(1 - 2en-V8-^)]2ri-V2(i + ^ (1 + ^ ^(^^ as n ^ 00. 

Finally, consider the contribution from terms with 121nn < A; < n — 5e(n + 3)^^^. By 
Lemma 2.4 this is at most 

9{n,k), 

12 In n<fc<n-5e(n+3)2/3 
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where g{n, k) = n^{n — k)\[{n + 2L)/e]^/n\. Now g{n, k+1)/ g{n, k) = {n + 2L)/{e{n — k)), 
which is increasing in k; hence g is U-shaped in k for fixed n (and L), with its minimum 
at 1 + |_n — (n + 2L) /ej . In particular, the maximum of g (n, k) for fixed n over an interval 
of values of k is achieved at an endpoint. Thus msix{g{n, k) : 12\nn < k < n — 5e(n + 
3)^/^} is achieved at an endpoint. At the first endpoint, k = [12 Inn] and g{n,k) = 



n e~ 



ri21nn' 



U.i=o[{n + 2L)/{n-i)] < n%-^^[{n + 2 L)/{n- 12 In n)]^^^^''''^ = (l + o(l))n- 



At the second endpoint, k 



n — 



g{n, k) 



n 



5e(n + 3)^/^ and the inequality n\ > (n/e)" gives 

[(n + 2L)/e]"-N"+3)^^'l/n! 
! [(n + 2L)/e]-r5-("+^)'''l (1 + 2Ln-i)"-N"+3)''1 
= (1 + o(l))e2^n6( \5e(n + 3)^/3] /e) l^'^i^+^m V27r5enV3[e/(n + 2L)] 1^^-+^)'^ 

^ [5e(n+3)V3] 



5e(n + 3)2/2' 
5e(n + 3)2/2' 



(l + o(l))e2V 



5e(n + 3)^ 



n + 2L 



V27r5en2/3 



by Stirling's formula. This clearly tends to faster than any power of n and is thus 
o(n~^). Now the terms at both ends have been shown to be at most (1 + o{l))n~^ 
and there are at most n terms, so the total contribution from these terms is at most 
(l + o(l))n-5 = o(l). 

To sum up, for n sufficiently large, X^fc f{n, k)[{n + 2L)/e]^/n\ is at most the sum of 
the contributions from the three ranges, which was computed to be (1 + o{l))e'^/{e — 
l)2 + o(l)+o(l). Thus by Lemma 2.1, P(0 is connected to i) > (1 + o(l))(e - 1)76^. □ 

Theorem 3.2 Let each edge of Bn be independently open with probability p = c/n. Then 
P(0 is connected to i by an oriented open path) converges to a limit as n ^ oo. The 
limit is if c < e and is (1 — x{c))'^ if c> e, where x{c) is the extinction probability for 
a Poisson(c) Galton-Watson process, namely, the solution in (0, 1) to x = e'^^^~^\ 



Note that as c ^ oo, x{c) = (1 + o(l))e 
probability is 1 — (1 + o(l))2e~'^ — >• 1. 



o(l), so that the limiting connection 



23 



Proof: There are n! oriented paths from 6 to i. Let N be the random number of these 
that consist entirely of open edges. For each path 7 the probabihty that 7 is open is p", 
so EA^ = n!p". If c < e then EN = nXd'n-'" = (1 + o(l))(c/e)''V27m, which tends to 0. 
Since P(A/" > 0) < EA'", this proves the first part. 

For the second part, fix c > e. Also fix e > 0. Write M — [1/e] . For %— 1, 2, . . ., let 
Ai be the set of vertices at level % reachable from in B„. For any fixed as n ^ 00, 
the joint distribution of l^ol? ■ ■ ■ 1 l^il approaches in total variation the distribution of a 
Galton- Watson process with the number of offspring of each particle Poisson distributed 
with mean c. Because a surviving branching process proliferates, an integer L = L(e) 
may be chosen so that P(|Ai| > M) > (1 — a;(c)) — e for sufficiently large n, where x{c) is 
the extinction probability for the Galton- Watson process, namely, the solution in (0, 1) 
to X = e'^^^~^\ Let Bj be Aj upside down, i.e., the set of vertices at distance j from i 
that can reach i. Then by symmetry and independence, we have P{F) > (1 — a;(c))^ — 2e, 
where F is the event {\Al\ > M and \Bl\ > M}. Now if either of the two sets A^ or 
is empty, then is not connected to i, so the limsup of (1 — a;(c))^ is established by the 
convergence in total variation. 

For the lower bound we employ the enhanced second moment method described fol- 
lowing Lemma 2.1, although the details here are slightly different. The variance-absorbing 
random variable Z is mind/lil, |-Bl|)- Uniformly in 2; > M, we show 

9 n-l 



P{N > 0\Z = z)> (1 -0(1)) 



1 + 



M(e-1)2_ 



(10) 



as n — > 00. Then 



and so 



P(A^>0) > (l-o(l)) 



1 + 



M(e-1)2 



P(F) 



liminf P(A^ > 0) > 



1 -1 



1 + 



M{e-iy 



[{l-x{c)r-2e]. 
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Letting e J, gives liminf„ P(A'" > 0) > (1 — a;(c))^, as desired. 

Henceforth tacitly conditioning on Z — z > M, we prove (10) by applying the 
second moment method to a truncation N' of N, as follows. First, reduce the prob- 
ability that any given edge falling between levels L and n — L is open from p — c/n 
to p = e/n. The obvious coupling argument shows that this diminishes N stochas- 
tically. Then let Xi, . . . ,xm be an enumeration of the first M vertices oi Al and let 
yi, . . . ,yM be an enumeration of the first M vertices of Bl, in some arbitrary order- 
ing of the vertices at levels L and n L, respectively. Let Ni be the number of open 
paths connecting Xi to yi, and consider A^' := Y^fLiNi. Essentially the same calcula- 
tions that showed EA^^/ (EA^)^ to be bounded in the proof of Lemma 3.1 will be used to 
show E ( AT') V (EAT') 2 < (l+o(l))(l + eV[M(e- 1)2]) via the two results 'ENf/iY.Nif < 
(l+o(l))eV(e-l)2, uniformly in and E(A/'iAr,)/[(EA/'i)(EAr,)] < (l+o(l)), uniformly in 
pairs i j. Indeed, these last two statements imply E{N'y = J2i EA^f + E(A^jAr,) < 
(1 + o(l))[eV(e - lymm^? + (l + o{l))E^M^m(EN,) < (l + o(l))M[eV(e - 
l)2](EAri)2 + (1 + o(l))M(M - l){ENiy < (1 + o(l)) (1 + eV[M(e - 1)^]) {EN'f. Now 
simply apply Lemma 2.1 to deduce (10) for A^' and hence for N. 

Working for convenience with p — e/{n + 2L) in Bn+2L rather than with p = e/n in 
Bn, and given distinct vertices xi and X2 in level L and distinct vertices yi and t/2 in level 

n + L, let Ni denote the number of open paths from Xi to yi, i = 1,2. We must show 
EA^iV(EA^i)2 < (l + o(l))eV(e-l)2 and E(A^iA/'2)/[(EAri)(EA/'2)] < l+o(l) as n ^ oo, 
uniformly in the choice of Xi,X2,yi,y2- Now the interval from Xi to yi is isomorphic to 
Bn, so A^^i has the same distribution as the total number of open paths, A^, in the proof of 
Lemma 3.1. Thus we immediately obtain EA^f/(EA^f) < {1 + o{l))e^ / {e — 1)^ (uniformly 
in Xi). For the other inequality, mimic the calculation from Lemma 3.1 to get 

E(^i^2) < j:uH{n,L,k,x„X2,yi,y2) (^)" < v F r„ r M + /„, 
(EA^i)(EA^2) \{^) n\\ k \ e J 

\ n+2L j 

Now break the sum into three pieces again, corresponding to values k < 12 Inn, k > 
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n — and all k in between. This time use Lemma 2.5 instead of Lemmas 2.2 - 

2.4, one difference being that the contribution for terms with < A; < 12 Inn is now 
o{{k + l)(n — k)\)[{n + 2L)/e]^ /nl, and thus the contribution from those terms and the 
unit contribution from the k = term sum to 1 + o(l). As before, the contribution 
from large k is o(l), and the contribution from the first term in the bound (3) for 
intermediate values of k is also o(l). We finish the proof of the theorem by showing that 
the contribution from the second term in (3) is also o(l): 

J2 [5e(n + 3)2/3j (2n^/8) N"+=^)'''l-i[(n + 2L)/e]'/n\ 

121nn<A;<n-5e(n+3)2/3 

= (l + o(l)) — =n^/^^ 



2^/2^ VnVs, 

which vanishes at a rate faster than any power of n. □ 



3.2 Oriented first-passage percolation 

Now consider oriented first-passage percolation (OFPP). Give each edge in Bn an up- 
ward orientation and assign independent, identically distributed random variables X^. 
with common density / to each edge e. The problem in OFPP is to determine the mini- 
mum value over oriented paths from to 1 of the sum along the path of the X^s. Under 
mild conditions on /, it turns out (Theorem 3.5) that this random minimum converges 
in probability to l//(0) as n — > oo. By multiplying every edge-passage time X^ by a 
constant, it can be assumed without loss of generality (provided < /(O) < oo) that 
/(O) = 1. 

As Aldous (1989) points out, use of the exponential distribution f{x) = e^^ simplifies 
some of the calculations involved but is not necessary. The following two lemmas, treating 
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the exponential and more general cases, respectively, produce large deviation estimates 
that correspond to the probability [e/(n + 2L)]^'*~^ in the proof of Theorem 3.2. 



Lemma 3.3 

(i) Let Sn be the sum of n > 1 independent, identically distributed random variables Yi, 
each exponential with mean 1, and let u be a real number in [0, 1]. Then P(<S'n < u) — 
(1 + Ki{u, n))e- V/n! with < Ki{u, n) < e/{n + 1) < 2. 

(ii) Given l<k<n-l, letS'^^ E^i Y, + Y^U+i Yl, where Fi, . . . , Y^, y.^i' ■ ■ ■ 
are independent and identically distributed. Then P(5'„ < 1 and -S^ < 1) < K2R{n,k), 
where 

Rin, k) = 22"-''=e'"-'=(2n - ky^^^'-''^ [{n - kf'''{2n - kf"' 

and K2 is constant. Furthermore, for 1 < k < n — 1, R{n, k — 1)/R{n, k) < K^/n for 
some constant K^. 

Proof: The key for (i) is the standard switching relation P(5'„ < u) = P(X„ > n), 
where X = (Xu)„>o is a Poisson process with unit intensity parameter. Thus P(S'„ < 
u) > P{Xu = n) = e^^M^/n!. Moreover, 



Using the switching relation together with this crude upper bound, we obtain 

P{Sn<u) = P(X„ = n)+P(X„ >n + l) 



n\ n-\-l' 
n\ n + 1 
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as desired. 



For (ii), we both prove the large deviations inequahty and show that it is tight. Begin 
by writing P{Sn, S'^ < 1) ^ /q P{Sk e du) \P{Sn-k < 1 - u)f. Using the bounds from 
(i) gives 

P{Sn,S' < 1) = /'[e-"H^-V(A;-l)!][e'"-2(l-H)'"-2*^][l+Xi(l-ii,n-A;)]V[(n-A;)!]2dii, 
Jo 

and using lower and upper bounds for e"~^ and Ki{l — u,n — k) bounds this below by 

{k-l)\[{n-kWJo ^ > 

and above by 

{k-l)\[{n~kWJo ^ > 

Now the integral is equal to [k — l)!(2n — 2k)\/{2n — k)\, so P(S'„, S'^ < 1) is bounded 
between and 9e~^ times (^^^~^'^^/(2n — /c)!. This can be approximated using Stirling's 
formula, for which it will suffice to note that the error factor of e^/^^^") is bounded. Thus 
P(5'„, 5"^ < 1) is bounded between positive constant multiples of R{n, k). 

To see that R{n,k — 1)/R{n,k) < K^/n, note that the exact quotient is 4e[(2n — 
k)/{2n -k + l)]2"-*^(2n -k + 1)-^ [{n -k)/{n-k + l)Y/^[{2n - k)/{2n -k + 1)]^/^. The 
factor [{n — k)/{n — k + 1)]^^^ is between and 1, while the rest of the product equals 

(1 + o(l)) 4/(2n — A;) < (1 + o(l)) 4/n, uniformly in k. This proves the claim and finishes 
that of the lemma. □ 

The Lipschitz condition in the following lemma does not give the most general / for 
which the first-passage times can be calculated, but it does cover most non-pathological 
cases. 

Lemma 3.4 Let f be a probability density on [0, oo), and suppose that /(O) = 1 and 
that f satisfies a "global" Lipschitz condition at the origin: \f{x) — /(0)| < K^x for 
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some positive K4 < 00 and all x > 0. Let Tn be the sum of n independent random 
variables Zi, . . . , Z„ with common density f and let Sn be the sum ofn i.i.d. exponentials 
Yi, . . . ,Yn with unit mean. Then for < u < 1, P(T„ < u) < e*^^"'"'^*'-*" P(S'„ < u) and 
liminf„P(T„ < u)/P{Sn < u) > e-^^\ Similarly, if T^^ = Eti + ELfc+i ^here 
Zi, . . . , Z„, . . . , are i.i.d., then P{Tn < 1 and < 1) < K^R{n,k), where 

Proof: For the upper bound, note that the Radon-Nikodym derivative of Zi with respect 
to Yi at X is f{x) /e-^ < {I+K4X) /e"^ < e^^+^^^''. Thus the Radon-Nikodym derivative of 
the n-tuple {Zi, . . . , Z^) with respect to (Yi, . . . , Y^) at (xi, . . . , Xn) is at most e(^+-'^4) E^^*, 
and hence the derivative of T„ with respect to Sn at x is at most e(^+^4)^. This estabhshes 
the upper bound for P(T„ < u). Together with (ii) of Lemma 3.3, this argument also 
estabhshes the upper bound for P(T'„ < 1,T^ < 1). 

For the lower bound, we first establish the fact that for any fixed u e (0, 1] and 
5 e {0,u), P(maxiy^ < S\Sn < u) converges to 1 as n — > 00. To see this, note that 
P{Sn < u,meiXiYi > 5) < nP{Sn < u,Yi > 5) < nP{Sn-i < u - 5); part (i) of the 
last lemma shows that for n > 2 this is at most 3ne^e~'^{u — 5Y~^/{n — 1)!. Then 
the lower bound from (i) of the previous lemma shows that P(maxi Y^ > 5 \ Sn < u) < 
2)e^{{u — S) / u)'"~^n^ / u which vanishes at an exponential rate as n ^ 00, proving the 
claim. 

Now, given K'^ > K4, a lower bound for the Radon-Nikodym derivative of Zi with 
respect to is 1 — K4X > e~^''i^l[x<s} for an appropriate 6 e (0, u). Then P(T„ < u) > 
e-K^p(^Sn < u,meiXiYi < 5); thus P(r„ < u)/P{Sn < u) is at least e--^>P(maXiyi < 
S\Sn < u) and so has a lim inf of at least e~^4" by the fact in the previous paragraph. 
Now let K'^ i K4. □ 

We are now ready for the main theorem for OFPF. 
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Theorem 3.5 Let the edges of Bn be assigned i.i.d. positive random passage times with 
common density f , and assume that — 1| < K4X for all x > 0. Then the oriented 
first-passage percolation time T — T*^"^ for Bn converges to 1 in probability as n —> 00. 

Proof: Let e be small and positive. With being i.i.d. with common density / and 7 
an oriented path from to 1 in B„, let Tn{'y) be the sum of X^y^ along edges vw of 7, so that 
the first-passage time T is just min-yT„(7). Let Tn denote a random variable distributed 
identically to each T'n(7) ^-nd let Sn denote the sum of n i.i.d. exponentials of unit mean, 
as in the lemmas. Showing that P(T < 1 — e) — > is easy. Let N be the number of 7 for 
which 7^(7) < 1 - e. Then, using Lemmas 3.4 and 3.3, P(A^ > 0) < EN ^ n!P(r„ < 
1-e) < n\exp[{l + K4:){l-e)]P{Sn < 1-e) < 3exp[X4(l-e)](l-e)'' = o(l) as n ^ 00, 
where Sn is the sum of n i.i.d. exponentials of mean 1. 

To show that P(T < 1 + e) is bounded away from 0, one can mimic the proof of 
Lemma 3.1, but in order to show that this probability converges to 1, we need to find 
another auxiliary random variable to reduce the variance. It will be easier to work 
in Bn+2- Let Aq be the random set of neighbors v oi for which the edge Ov has 
< e/2. Similarly, let Ai be the random set of neighbors v oi 1 for which the edge 
vl has X^i < e/2. Let b' be the minimum of \Ao\ and \Ai\. Enumerate the elements of 
Aq by Xi,X2, ■ ■ ■ and the elements of Ai by yi, y2, ■ ■ ■ in such a way that for 1 < i < b', 
yi lies above Xi. This is easy to do since there is only one neighbor of 1 that does not 
lie above any given Xi. Let b — [en/4]. The first thing to observe is that b' > b with 
probability converging to 1 as n ^ 00. This is immediate from the fact that b' is the 
minimum of two independent random variables that are binomial with parameters n-\-2 
and Jo^^ f{x) dx > Jq"^ (1 — K^x) dx > e/4. Now condition on the event that b' > b. It 
suffices to show that the probability of finding an oriented path 7 connecting Xi to yi 
with T„(7) < 1 for some i < b converges (with appropriate uniformity) to 1, conditionally 
given b' > b and the enumeration of the Xj's and t/j's. What will in fact be shown is 
that, uniformly over all choices of vertices xi, X2, ■ ■ ■ ,Xb neighboring 6 and yi, 7/2, • • ■ , Z/6 
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neighboring i with yi above Xi for each i, the probabihty of finding an oriented path of 
passage time at most 1 connecting some Xi to yi tends to 1. 

For this we use the second moment method. For 1 < i < 6, let be the number of 
paths connecting Xi to yi with passage time at most 1. Let N — Yh=i ^i- The interval 
in Bn+2 from Xi to yi is isomorphic to B„. It is therefore easy, using Lemmas 3.3 and 3.4, 
to see that 

b 

EN = ^EiVi = 6EiVi = 6n!P(T„ < 1) = bco, 

i=l 

where Cq = Co(n) is bounded between positive constants. Now 'EN'^ — X^^EA^? + 
E(NiNj). If we can show that 

EN^ = 0(1) and E{N^N2) < (1 + o(l))c^, (11) 

uniformly in the choice of xi, X2, yi, y2, then, it will follow that EA'"^/(EA'")^ < 0{b~^) + 
1 + o(l), uniformly in the choice of the 2b vertices. Since b tends to infinity with n, this 
bound converges to 1, and so F{T < 1 + e) — > 1, proving the theorem. 

Each part of (11) is established in pieces, in a manner similar to the bounding of (9). 
For any fixed 7 connecting xi to yi, ENf is given by nl times the sum over 7' connecting 
Xi to yi of P{Tn{'~f) < 1 and T„(7') < 1). Break the sum into three ranges according 
to the number k of edges shared by 7 and 7' as before, and additionally separate the 
cases A; = and k — n. The case k — means Tn{'y) is independent of 2^1(7'), so the 
contribution to EA^^in this case is at most (EA^i)^; and the case k — n has 7 = 7', so 
the contribution in this case is exactly EA^i. Using Lemma 3.4 for the other three ranges 
and recalling that the case A; = n — 1 is impossible, the sum can be bounded by 

n-2 

ENf < {ENif + ENi + n\^f{n,k)K5R{n,k) 

k=l 

< (co + co) 

+ (l + o(l))n! {k + l){n-k)\K^R{n,k) 

0<fe<121nn 
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+ n! Yl (n-k + l) {2n"^Y-^K^R{n, k) 
+ n! n^{n-k)\K^R{n,k), 

12 lnn<A;<n-5e(n+3)2/3 

for large enough n. Now it will not be too hard to show that this is Cg + cq + 0(1) = 
0(1), but before doing so, notice how similar the above bound is to a good bound on 
'Ei{NiN2). The A; = term for Yi[NiN2) has the same bound as above, but the k — n 
term vanishes. The necessary changes are completed by using Fi in place of /. This 
allows the {k + l){n — k)\ to be replaced by o{{k + l)(n — k)\) according to Lemma 2.5, 
and hence 



E(A^iA^2) < cl + oln\ Y {k + l){n-k)\K^R{n,k) 

\ 0<fc<121nn 

+ n! {n-k+l) {2r{"y-^K^R{n, k) 

n-2>A;>n-n3/4/2 

+ n! Y 2n\n-k)\K^R(n,k) 

12 In n<fe<n-5e(n+3)2/3 

+ n! Y f5e(n + 3)2/3] (2^^/^) N"+3)'/'l-iK5i?(n, /t). 

12 In n<fe<n-5e(n+3)2/3 

It will be shown that the last three terms of the bound on EA^f are respectively 0(1), 
o(l), and o(l). This will show that EA'"f = 0(1), and also that the first four terms of the 
bound on E(A''iiV2) sum to Cq + o(l) + o(l) + o(l) = Cg + o(l). For large enough n, the 
fifth term of the bound on E(iViA^2) is bounded by 



;i + o(l))i^5('^V^'n^^/^^ 



X 1 + 



5e(n + 3)2/3 



(n+|"5e(n+3)2/3]) 



8e \ r5^("+3)^/^l 



n 



which vanishes at a rate faster than any power of n. Thus E{NiN2) = Cq + o(l), com- 
pleting the proof of the theorem via (11). 
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The three estimates for the bound on ENf are now routine calculations. Plugging in 
the value of R{n, k) and using Stirling's formula gives for the second term in the bound 

(l + o(l))n! ^ {k + l){n-k)\Kr,R{n,k) 

0<fe<121nn 



l<fe<121nn 



= (1 + o(l))X52'/'7r + [(n-A;)"-V22"-2fe(2n-A;)-(") 

l<fc<121nn 

Note that if the sum here were to contain a A; = term, that term would equal 1. 
Furthermore, changing A; to A; + 1 multiplies the part of the summand inside square 
brackets by 

{n-k- l)-^[{n -k-l)/{n- A;)]"-*^2-'[(2n -k- l)/{2n - A;)]-('"-'=)(2n -k-l). 

Now [{n-k - l)/(n - k)f-^ e-^ while [(2n - k - l)/(2n - A;)]-^^"-*^) e and 
{2n — A: — l)/(n — A; — 1) — > 2, all uniformly over k < 12 Inn. Thus the successive ratios 
arc (1 + o(l))| uniformly over k in the range of summation. Therefore the sum is at most 
(1 + o(l)) E^i(^ + IjS"'^ = (1 + o(l))3 = 0(1), establishing the first bound. 

For the third term in the bound on ENf, let m = n — k. Then plugging in for R[n, k) 
and using Stirling's formula yields 

n\ Yl (n-k + l) {2n''I^Y-^K^R{n, k) 

n-2>fe>n-n3/4/2 

= {l + o{l))K^ Y (n- A; + l)(2n^/^)"-*^n"e-"v^22"-2V"-'= 

n-n3/4/2<A;<n-2 

x{2n - A;)-(^"-'=)(n - A;)-^/^(2n - k)-^/^ 

2<m<n3/4/2 
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< (l + o(l))i^5(27r)^/2 ^ 2m(8e)™r^-"^/^ 

2<m<n3/4/2 

For the first inequality here we used ^Jn(^n — k)^^!"^ = 1 + o(l) uniformly over k > 
n — and for the second we used m^/^ + m~^/^ < 2m and (n + m)~^"+'"'* < 77,-("+'"). 

Changing m to m + 1 multiplies the term by 8e(l + l/m)n~^^^, which vanishes in the 
limit uniformly in m; thus the sum is dominated by the m = 2 term, whose value is a 
constant times and is thus = o(l). 

Finally, to bound the fourth term in the bound on ENf, plug in to get 

n\ Yl n%n-k)\K5R(n,k) 

12 lnn<fc<n-5e(n+3)2/3 

= (l + o(l))27rii'5 Yl n^(n- A;)"-*^e-("-^')y^^r^n"e-"v^22"-2'=e" 

12 In n<fe<n-5e(n+3)2/3 

x(2n - k)-^''^-'\n - k)-'/\2n - A;)-^^ 
< (1 + o(l))27rii:5 Y n^(n-A;)"-*^n"2'"-''=(2n-A;)-("), 

12 In n<fe<n-5e(n+3)2/3 

since ^/n{2n - k)~'^/^ < 1. 

The sum here is at most n times its largest term. Let r = k/n and rewrite the typical 
summand as n^[(4 — 4r)^~''/(2 — r)2-rjn _ j^e^j-^^^n^ gg^y_ Now we find the maximum of 
h{r) on [0, 1]. Taking logs gives 

ln/i(r) = (1 - r) ln(4 - 4r) - (2 - r) ln(2 - r), 

so that 

(d/dr) In h(r) = ln(2 - r) - ln(4 - 4r). 

This increases from — ln2 to oo as r increases from to 1, so the maximum of \nh{r) 
over the interval {12lnn)/n < r < 1 — 5e(n + 3)^/^/n is achieved at one of the endpoints, 
at least for large n. Again for large enough n, we can, for any S > 0, get the derivative 
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of Inh on [0, (121nn)/n] to be bounded above by — ln2 + 5 < —0.693 + S, so choosing 
S = 0.003 makes the derivative of In h bounded above by —0.69 on this interval. Similarly, 
for large enough n the derivative on [1 — 5e(n + 3)^/^/n, 1] is bounded below by 1. 
Noting that h{0) = h{l) = 1, it follows that the value of n^h{r)^ at r = (12 In n)/n 
is at most [e(-0-69)(i2inn)/njn^6 ^ ^"2.28^ ^^^^^ at 1 - 5e(n + 2>f/^/n is at most 

j-g-5e(n+3)2/3/njn^6 _ g-5e(n+3)2/3^6 ^ j7,-2-28 large u. Thus the sum Under consideration 
is at most 

Putting all of this together gives EiV^^ < cg + co + 0(l) + 0(n-^/^) + 0(n-^-28) = 0(1), 
as desired. □ 



4 Unoriented percolation 

In Section 5 we shall consider the first-passage time to 1 for unoriented first-passage per- 
colation on Bn- For completeness, in this section we treat ordinary unoriented percolation 
and argue that the critical probability is as put forth in the following theorem: 

Theorem 4.1 Let each edge of Bn he independently open with probability p = c/n, < 
c < oo. Then P(0 is connected to 1 by an (unoriented) open path)^ (1 — x{c)Y , where 
x{c) is, as in Theorem 3.2, the eoctinction probability for a Poisson(c) Galton-Watson 

process. 

Proof: Write 9.^ = 0n{c) for the percolation probability in question. We first note that 
limsup„^„ < (1 — a;(c))^ by a branching process approximation similar to that in the 
second paragraph of the proof of Theorem 3.2; we omit the details. 

For the lower bound we may restrict attention to the case c > 1; it is precisely for 
these values of c that y{c) := 1 — x{c) > 0. Let < e < y{c)/A. We rely heavily on a 
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result of Ajtai, Komlos, and Szemeredi (1982): P{F) > 1 — o(l), where F is the event 
that (a) there is exactly one component in the random graph formed by the open edges 
that has at least (y(c) — e)2" vertices, and (b) all the other components are of size at 
most e2". 

For V & Bn, let Ay denote the event {y is connected by an open path to at least 
{y{c) — 2e)2" vertices}. By a simple application of the FKG inequality (Fortuin, Ginibre, 
and Kasteleyn (1971)), the indicators of the events arc pairwise positively correlated. 

Furthermore, conditionally given F, we have by symmetry 

P{Ay) = P{v e the unique giant component (GC)) = 2""E(size of GC) > y{c) - e; 
thus, unconditionally, P(^^) > y{c) — 2e for sufficiently large n. By FKG, 

P(^6n^i)>(z/(c)-26)^ 

and so 

P(0 and i are in the same component | -F) > P(6 and 1 are in the GC | F) 

. P(^, n A, I F) > P(-^.nA.)-(i-P(f)) ^ j^^^^ _ 3^j,_ 

and hence 6n > {y{c) — 4e)^, for sufficiently large n. Let e J, to complete the proof. □ 

We close this section by noting that for c < 1 there is a more elementary proof that 
^n(c) — > 0. For X e Bn, let g{x) denote the probabihty that x is connected to by an 
open path. Clearly, for x ^ 0, g{x) equals the probability that there is a neighbor y of 
X such that y is connected to by an open path not containing x and the edge {y, x} is 
open. Hence 

5(0) = 1> 9i.x) < PJ2y^x9{y) for Xy^O, g{x) < 1 for all x, (12) 
where the sum is over vertices y adjacent to x. 

Repeatedly applying (12), we find g{x) < pn = c for x 6, g{x) < for (i(0, x) > 2, 
g{x) < for d{Q, x) > 3, . . . , and finally 6n = g{l) < c" — >• 0, as desired. 
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5 Richardson's growth model and unoriented perco- 
lation 



Consider the following model for the spread of disease. Individuals are located at ver- 
tices of an n-cube, with edges modelling pairs of individuals in frequent contact. One 
individual, 0, is infected at time and the rest are healthy. Independently for each edge 
between an infected individual and an uninfected one, there is a constant small proba- 
bility per small unit of time that the contact between those two individuals will cause 
the uninfected one to become infected. It is easy to construct from this description a 
stochastic model for the growing set of infected individuals. The model is a continuous 
time Markov chain on the space of subsets of Bn which jumps from A to AU {v} at 
rate k{A,v), where k{A,v) is the number of infected neighbors of v, i.e., the number of 
neighbors of v in A. This Markov chain is called Richardson's growth model. 

Interesting questions about this model are (1) When should we expect 1 to become 
infected? and (2) What is the cover time, i.e., when should we expect all the vertices 
to become infected? In this section we discuss the first question, giving limiting upper 
and lower of 1 and 0.88, respectively. These are obtained by proving and then exploiting 
the fact that the infection time for i in Richardson's model has the same distribution as 
the first-passage time to i in unoriented first-passage percolation on The cover time 
question is addressed in Section 6. 

The following lemma reduces the problem of when 1 first becomes infected to unori- 
ented first-passage percolation with exponentially distributed edge-passage times. Since 
the oriented percolation time is always at least as great as the unoriented percolation 
time (the minimum over paths directed away from must be at least the minimum over 
all paths), the upper bound of Theorem 5.2 for the infection time of i is immediate. 

Lemma 5.1 Let the edges {v,w} of the undirected graph Bn be assigned independent 
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exponential random variables X^^^ = Xw,v of mean 1. Define the infection time of a 
vertex v, denoted Tniy), to be mi^^X^^^^^^, where the inf is over all paths from f o = 6 
to V. Let A(t) — {v & Bn : Tn{v) < t}. Then the random map A from [0, oo) to subsets 
of Bn has the same law as Richardson's model. 

Proof: See Durrett (1988, page 177) for a sketch of this proof. □ 

Theorem 5.2 For any e > 0, the probability of finding 1 infected by time 1 + e in 
Richardson's model on Bn, beginning with only infected at time 0, tends to 1 as n ^ oo. 

Proof: Theorem 3.5 and Lemma 5.1. □ 

We doubt whether this result is sharp, since there is no reason why the unoriented 
percolation time should be as great as the oriented percolation time. The next theorem, 
based on a calculation by R. Durrett (personal communication), gets a lower bound 
for the unoriented first-passage time by comparing to a branching translation process 
(BTP). This is a process, started with a single particle at 6, for which each existing 
particle generates offspring at rate n, where the offspring are each displaced from the 
parent by an independent uniform random step e(j). Letting Z{x,t) be the number of 
particles at x at time t, the process {Z{x,t) : x G K„, t > 0) is formally defined by 
the transition rates {Z{x)) —>■ {Z{x) + 6xy) at rate J2w:d{y.w)=i Z{w), where 6xy is 1 if 
X = y and otherwise. It is easy to couple BTP to Richardson's model so that the set of 
infected vertices in Richardson's model is always a subset of the set of populated vertices 
in BTP. Thus the first time r„ that i is populated in BTP is stochastically less than the 
first infection time T„ of i in Richardson's model. 

Theorem 5.3 (Durrett) As n ^ oo, the time Tn of first population ofl in BTP con- 
verges m probability to ln(l + v^) = 0.88. Consequently, P(T„ < ln(l + v^) - e) ^ 0. 



38 



Proof: The lower bound will be gotten by a routine first moment calculation. The upper 
bound in probability for BTP (which is not necessary for the result on Richardson's 
model) requires a second moment calculation and a little more work. Fix n and write 
mi{x,t) for 'EiZ[x,t) = the expected number of particles at x at time t in BTP starting 
from a single particle at 0. We remark for later that this is also the expected number of 
offspring at y + x at time s + t oi a particle at y at time s that are born to the particle 
after time s, where the addition in y + x is taken, as usual, to be coordinatewise mod 2 
addition. Since P(Z(1, t) > 0) < mi(l, t), the lower bound in probability will follow from 
showing that mi(l, — > as n — > oo for any i < ln(l + ^/2). Viewing vertices of B„ as 
sets, we write \x\ for the cardinality of x; the differential equation for mi(x, is easily 
seen to be 

— — = 2^ m{y,t) (13) 

y:d{x,y)=l 

with initial conditions 'm{x, 0) = S^^. Let 



p{x,t) = 




be the probability that a simple random walk with rate n started at is at x at time t. 
Then, as may be verified by a variety of probabilistic and analytic arguments, the unique 
solution to (13) is given by 

mi{x,t) = e^^p{x,t). 

Putting X = i gives 



mi{i,t) = 




Since (e* — e *)/2 is increasing in t and equal to 1 at t = ln(l + \/2), it follows that for 
t < ln(l + \/2), mi(i, t) tends to as n — >• cxd. Hence P(r„ < t) — >■ 0, as desired. 

The upper bound in probability on t„ is gotten by a now familiar sort of argument. 
Fix e > 0. First the second moment method is used to show that liminf„P(T„ < 
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ln(l + v^) + e) > 1/12. Then the initial branching of the process is used to show that 
with just 2e more time units, there are actually many independent chances of no less 
than 1/12 each to get 1 populated, and hence the probability that this occurs is near 1. 

Begin with 

(EZ(i,t))2 mi(i,t)2 



P(T„<i) =P(Z(l,i) >0) > 



E(Z(l,t)2) m2(l,t) 
where m2{x,t) := E{Z{x,t)'^). To calculate the value of m2{x,t) in terms of mi{x,t), 
write Z{x,ty as Z{x,t) plus twice the number of unordered pairs of distinct particles at 
X at time t. Each such pair of particles has a well-defined time s at which their ancestral 
lines first split apart. Say that at time s a particle pi at vertex y gave birth to a particle 
P2 at vertex y + e{i), and that both particles are descendants of pi but only one is a 
descendant of p2. For fixed y and i and interval [s, s + ds), the expected number of such 
pairs is mi{y, s) dsmi{x — y,t — s) mi{x — y — e{i),t — s), so summing over y and i and 
integrating over s gives 

1712(1, t) 1 ft mi{y,s)mi{i-y,t-s)mi{i-y-e{i),t-s) 



mi(l,t)2 mi(l,i) i Jo y mi(l,i)2 

(14) 

Now fix t = ln(l + ^/2) + e. The first term tends to as n — >• oo, so it suffices to show 
that the limsup of the sum on i is at most 12. 

Substituting ■mi{x,t) = e"'^p{x,t) into the sum on i in (14) yields 

2 r do y c""' ^^^' ^^^^^ -y^'t- s)p{i -y- e{i),t - s) 

Next, plug in the value for p{x,t). At the same time, group together all y on the same 
level of Bn, i.e., all y with \y\ = k for each k. Then \y + e(i)| will equal either A; + 1 or 
k — 1; since p{l — x,t) increases with we get an upper bound by replacing \y + e{i)\ 
by A; + 1. This gives an upper bound for the integrand of 

ne — 



-2{t-s) 
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X E 2-" r [(1 - e-2^)(l + e-2(*-^))2]'=[(l + - e-2(*-))2]-*^(i _ e-^*)-^- 

fc=o V^y 

The sum over k is just the binomial expansion of 

/ (I - e-^^)(l + e-^(*-^))^ + (1 + e-^^)(l - e-^^*-^))^ " 
V 2(1 - e-2*)2 J ' 

and simphfying this yields 

which gives a bound for the integrand in (15) of 



1 + e-'2{t-s) ^ / g-4(t-5) 



e 



-4t ^ 



1 + e"2(*"^) 

We need a better bound on the integrand when s is near t: the factor — — r 

blows up like {t — s)~^, which is not integrable. Note that in the case A; = n it is not 

possible to have \y + e{i)\ = k + 1. Thus for the k = n term, the factor _2{t~s) 

be replaced by its reciprocal. This reduces the integrand significantly when the k = n 
term is the dominant term in the sum. The ratio of the k = n term to the entire above 
sum on k is 

[(1 - e-'^)(l + e-2(*-^))2]'*/[(l - e-'^)(l + e-^^'-'^f + (1 + e-^'){l - e-^^*-^))^]" 

> [(1 - e-^^)/((l - + (1 + e-^'){2{t - s))^)]" 

> [l-X(t-s)2]" 

> l-K{t-s) 



for some constant K when t — s < 1/n. {K — 15 will do when n > 2.) Also, 

^1 _ e-2(t-s)- 



[ _o(t_ J ^ ^ {t- sY + 0{{t - s)^). Putting this all together, a better bound for 
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the integrand in (15), uniformly for s satisfying t — s < 1/n, is 



[l + o{l))K{t-s)n 



1 + e-^(*-^) 

1 _ p-2{t-s) ' 



1 + 



g-4(t-s) _ g-4t ^ 



(1 



(17) 



In order to make use of (16) and (17), examine the function 

g-4(u-s) _ g-4u' 



G{s, u) = In 



1 + 



(1 - e-2")2 

This is convex in s for the values of u we are interested in, which may be seen by 
differentiating twice with respect to s: writing c = c„ = e~^"/(l — e"^")'^ gives 

d^G{s,u) 



d 



-s + ln(l + c(e^^-l)) 



l + c(e4* - 1)^ 
16c(l - c)e^' 



[l + c(e4--l)]2' 

which is positive for all s whenever c < 1. This is true if and only if m > In a/2 and hence 
for u — ln(l + ^/2) + e. In particular, the maximum of G(s, u) over s £ [0, u\ is achieved 
at an endpoint. But G(0, u) — and 

1 - e-^" 



G{u, u) = In 



6-^1 + 



In 



2e"/(e^" - 1) 



(l-e-2")2' 

This is decreasing in u and has value when u = ln(l + \/2). Thus when u = t = 
ln(l + v^) +e, G{u,u) has a negative value which we shall call —V{e). Now we bound 
the upper bound for (777,2(1, t) — 777i(i, t))/777i(i, t)^ given by display (15) in three pieces: 



< 2 



7772(i,t) - 777i(i,i)]/r77i(i,0' 

1/2 I _^ g-2(t-s) 



1 - e-2(*-*) 



(18) 
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A/2 1 - e-2(*-^) ^ ' 

+ (1 + o(l))2 - (20) 

For the first piece we calculate the value of G'(l/2,t), getting a constant less than 

-1/4. Thus by convexity, G{s,f) < -s/2 for < s < 1/2. Now for < s < 1/2 

1 + e-2(t-6-) 

and any e, the factor -77 — r is at most 3, so the contribution from (18) is at most 

\ — g 2(t s) 

2/0/^ 3ne-"^/2ds = 12(1 - e""/^) < 12. 

For the second piece, bound G{s, i) on 1/2 < s < i by its values at the endpoints; for 
small e the greater value is the value at the right endpoint, namely, —V{e). The value 

of -77 — r onl/2<s<t — 1/nis bounded by its maximum, which is achieved at 

\ — g 2(t s) 

s — t — 1/n and has a value of at most n + 1. Thus the contribution from (19) is at most 
2 /i/V^" + l)^e-"^(^) ds < 2t{n + l)ne-"^(^), and this tends to as n ^ 00. 

To bound the third piece, expand the integrand in powers of (t — s) to compute the 
integral as (l + o(l)) {t- s)-^nK{t- s)e''^^''^^ ds < {l + o{l))Ke-''^^^\ which tends 
to as n — > 00. Thus the entire integral is bounded by 12(1 + o(l)), and liminf„ P(t„ < 
ln(l + V2) + e) > 1/12. 

Finally, to show that P(r„ < ln(l + a/2) + 3e) 1, let A be the set of particles at 
distance 1 from at time e. Then \A\ > n^/^ with probability approaching 1 as n ^ 00. 
The offspring of elements of A now act independently from time e to time ln(l + \/2) + 2e, 
each particle at some y with d{6, y) = 1 having probability at least (1 + o(l))^ of having 
a descendant at the antipodal point to y at time ln(l + \/2) + 2e, according to the 
calculation just completed. Letting B be the set of particles at sites that neighbor i at 
time ln(l + v^) + 2e, P{\B\ > n^/^) is at least P{\A\ > n^/'^)P{X > n^/^), where X is 
a binomial with parameters n^/^ and (1 + o(l))^. This probability also tends to 1 as 
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n-^ oo. Finally, P(r„ < ln(l + V2) + 3e \ \B\ > n^/^) > 1 - e"''*'^", which tends to 1 as 
n ^ oo, proving the theorem. □ 

6 Covering times in Richardson's model 

This section answers affirmatively the question of whether the time until the entire n- 
cube Bn is infected is bounded in probability as n — > oo. The constant upper bound 
given here is 41n(4 + 2-\/3) + 6 = 14.04. This is by no means sharp, but on the other 
hand we also produce a lower bound on the cover time of | ln(2 + + ln2 = 1.41. 
Since the infection time of 1 is bounded between ln(l + -\/2) = 0.88 and 1 in probabihty, 
this means that the liminf in probability of the cover time is strictly greater than the 
limsup in probability of the time to reach the farthest vertex. Perhaps the cover time 
has a limit in probability, but we do not venture a guess as to what the limit should be. 

6.1 The upper bound 

The following statement of duality in Richardson's model will be helpful. The proof 
can be found in any introduction to the contact process, such as Durrett (1988). The 
intuition is to think of A2{t — s) as the set of vertices that would be able to infect 1 by 
time t if they were infected at time s. 

Lemma 6.1 Let {Ai{t)) be Richardson's model on Bn as defined above and let {A2{t)) be 
an independent copy with the difference that ^2(6) is set to be {1} instead of {0}. Then 
P(i e Ai{t)) = P{Ai{s) n A2{t - s) 7^ 0) for anyO<s<t. □ 
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Theorem 6.2 Let {A(t)) be Richardson's model on Bn starting with only infected at 
time 0. Let c := 21n(4 + 2y/3) + 3 = 7.02. Then for any e > there exists N = N{e) 
and h{e) > such that for n>N, P(i e A{c + e)) > 1 - (4 + h{e))-''. 

We defer the proof of Theorem 6.2 in order to present the resulting cover time upper 
bound. 

Corollary 6.3 For any e > 0, P(A(2c + e) = Bn) 1 as oo. 

Proof of Corollary 6.3: Let = N{e/2) be chosen as in the statement of the previous 
theorem and pick any n > 2N. Let x be any element of Bn- First suppose that the 
distance from 0toa;is(i>n/2. Then the sublattice with top element x and bottom ele- 
ment is a Boolean algebra of rank at least N, and the induced process on the sublattice 
(i.e., the process with no infections allowed except on the sublattice) is still Richardson's 
model. By the previous theorem, x is infected by time c + e/2 with probability at least 
1 - (4 + /i(e/2))-<^ > 1 - (4 + h{e/2)y/'^. 

On the other hand, suppose x is at distance less than n/2 to 6. Then after time 
c + e/2 the top element i is infected with probability at least 1 — (4 + h(e/2))~'"' and, 
conditioned on that, the probability that x is infected another c + e/2 time units later 
is at least 1 — (4 + /i(e/2))~"'/^ by the previous argument (since the distance from x 
to i is at least n/2). Thus each x fails to be infected at time 2c + e with probability 
at most (4 + /i(e/2))^"' + (4 + /?.(e/2))~"/^, and summing over all x gives at most (1 + 
/i(e/2)/4)-"/2 + 2"", which tends to as n ^ oo. □ 

Notice that the reason we get 2c instead of c as an upper bound in probability for the 

cover time is that the proof of Corollary 6.3 gives better upper bounds on the probability 
that a vertex is uninfected the further it is from 1. It is unlikely that the bounds 
reflect the true state of affairs. In particular, we suspect that the random time until 
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a vertex x is infected is stochastically increasing in |a;|. This would immediately imply 
an upper bound in probability of c = 7.02 for the covering time. (In fact, the bound 
could then be lowered to c := 21n(2 + -\/2) + 3 = 5.46 by establishing the modification 
P(i e A{c-\-e)) > 1 — (2 + ^(e))~" of Theorem 6.2.) More generally, we have the following 
conjecture. 

Conjecture 1 Let G he a graph with distinguished vertex x. Let H he the graph G x 
{0, 1}; with edges hetween {x,i) and {y,i) for neighhors x,y of G and i — 0,1 and edges 
hetween {y, 0) and {y, 1) for all y & G. For z & H, let T{z) he the time that z is first 
infected in Richardson's model on H heginning with a single infection at (x, 0) . Then, for 
any y & G, T{y, 0) is stochastically smaller than T{y, 1). 

We now turn to the proof of Theorem 6.2. 

Proof of Theorem 6.2: Throughout the proof we use the notation x' for the complement 
of a vertex x (viewing x as a subset of {!,..., n}) and S' for {x' E Bn ■ x E S} 
when S C B„. Let to = 0, = ln(4 + 2^) + e/7, t2 = ln(4 + 2^3) + 2e/7, ^3 = 
ln(4 + 2^3) + 1 + 3e/7, and ^4 = ln(4 + 2V3) + 2 + 4e/7. By Lemma 6.1 it suffices to 
find h for which P{Ai{t4) n ^12(^3) 7^ 0) > 1 - (4 + /i(e))"" for large n. The method will 
be to watch the evolutions of Ai and A2 and look for vertices x for which simultaneously 
X e Ai{t) and x' e A2{t). In particular, we will show that for i — 1,2, 3, certain subsets 
Di oi {x : X E Ai{ti) and x' E ^2(^1)} (actually, of a slight modification of this for D3) 
are sufficiently large, and then we will argue that each x in D3 has an independent chance 
of becoming an element of ^1(^4) fl ^2(^3)- 

For 1 < i,j < n let rjij — {i, j} E Bn- Let rii — [n — 2 log2 n\ and let 5*2 C Bn be the 
set {rjij : i < j < rii}. Let Si — {{i} : i < rii} be the set of elements at level 1 of B„ 
beneath 5*2. For each x E S2 let T^. C Bn be the set { y : xn{l, . . . , rii} — yn{l, . . . , rii} }. 
Note that and Ty are disjoint for distinct x,y E S2; similarly for and Ty. 
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Several of the arguments below will involve the monotonicity of Richardson's model: 
forbidding some edges to pass the infection at various times decreases A{t) and hence can 
only increase the infection time to any vertex. It can therefore only increase all infection 
times to suppose for Ai that from time to to time ti infections may occur only in rank 1 
of Bn, from time ti to time t2 infections may occur only in rank 2 (and in fact later the 
set of allowed infections will be further restricted) , and from time ^2 to time ^3 infections 
may occur only between x and y when the symmetric difference x Ay is the singleton 
{i} for some i > rii (in other words, infections from t2 to may occur only when x 
and y are neighbors both in for some z). Also suppose dually for A2 that infections 
between times to and ti occur only in level n — 1, that infections between times ti and t2 
occur only in level n — 2, and that infections between times ^2 and occur only between 
elements of the same T^. Finally, suppose that between times and t^, infections in Ai 
spread only between vertices x and y for which x fl {rii + 1, . . . ,n} — yH {rii + 1, . . . , n}. 

Let Ui = SiC] Ai{ti), let Vi = S[n A2{ti), and let Di = UiH ¥[. The first claim is 
that there exist Si{e) and hi{e), both positive, for which P(|Di| < 6i{e)n) < (4 + /ii(e))~"' 
for sufficiently large n. This is just a large deviation calculation. \Di\ is the sum of ni 
i.i.d. Bernoulli random variables, each equalling 1 with probability p = (1 — e~*i)^ and 
with probabihty q — 2e~*i — e~^*i . Let < a < p and write b — 1 — a. By choosing the 
optimal value In(^) for 9 > 0, we find from the moment generating function inequality 

P(|Di| < am) = P(e-^l^il > e-^'^"!) < e^""iEe-^l^il 

that P{\Di\ < ani) < [{p/aY{q/bff\ As a i 0, {p/ayiq/bf I q = 26'^' - e^^*! < 
2(4 + 2 V3)-i - (4 + 2V3)-2 = (2 - v^) - (7/4 - ^3) = 1/4. Thus there exists 6o{e) > 
and /io(e) > such that P(|-Di| < So{e)ni) < (4 + /io(e))~"^ Since ni/n 1, this implies 
the existence of /ii(e) such that for any fixed 5i < 5o{e), P(|-Di| < 6in) < (4 + /ii(e))~"' 
for sufficiently large n. 

Now let U2 = 820 Ai{t2), let V2 = S'^ n A2(t2), and let D2 = U2 n The next 
claim is that there exist S2{e) and /12(e), both positive, for which P(|£^2| < S2n'^) < 
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(4 + /i2(e))~" for sufficiently large n. Between times ti and t2 the only infections allowed 
in Ai involve vertices in 5*1 infecting vertices in 5*2. For convenience, restrict further the 
allowed infections by requiring that {i} may infect {i.j} only iii < j and j — i < ni/2 or 
i > j and i —j > ni/2. Note that each {i,j} in can be infected by i or j but not both. 
The exception is \ j — i\ = ni/2; in that case {i.j} cannot be infected at all. Then for each 
X & Di, there is a set S{x) of [rii/2] — 1 vertices in 5*2 that x can infect between times 
^1 and t2, and these sets are disjoint as x varies over Si. For each x & Di, the number 
of vertices infected by x by time t2 whose complements have been infected by x' by 
time t2 in the process A2 (with the dual restrictions) is a binomial random variable with 
parameters [ni/2] — 1 and (1 — e~^/^)^. Now the probabihty that this binomial is less than 
half its mean is exponentially small in n, say < e""", so conditioning on \Di\ > 5in, the 
probability that no more than \5iri \ /2 of these i.i.d. binomials are greater than half their 
means is at most 2r''i"le~"'^i"^/^. This is smaller than (4 + /ii(e))~" for sufficiently large n, 
and thus we have shown that P (1^2 1 < {l-e-'/^f \Siri\ (K/2] - 1) /4) < 2(4+/ii(e))-" 
for sufficiently large n. Choosing S2 < (1 — e~^/^)^5i/8 and /i2 < hi proves the second 
claim. 

Condition until the last sentence of this paragraph on D2. Between times t2 and ^3 the 
spread of infection in A 1 is confined to each T^, so the spread of infection is independent on 
each Tz. The same goes for the propagation of A2 on each T^. On each T^, the process is 
just a Richardson's model on a cube of dimension n — rii; hence uniformly for x & D2 and 
y with ni) < d{x, y) < rii), the probability that y is infected by time ^3 = 
^2+(l+e/7) tends to 1 for large n by Theorem 5.2. Similarly, the probability that y is dual 
infected by time ^3 tends uniformly to 1, where y is the element of T'^ that agrees with y in 
the last n — ni places. (Note d{x' , y) > |(n — ni).) Thus the probabihty that y is infected 
by ^3 and y is dual infected by ^3 also tends uniformly to 1. Thus uniformly for x & D2, the 
expected cardinality oi Vx :— {y E T^; : y is 74i-infected by ts and y is 742-infected by t^} 
is at least (1 — o{l))\Tx\ > (1 — o{l))n^. Then, since we always have < \Tx\, P(|V^a;| < 
nV2) < 1/2 for large n, so P(|V;| < for all x e D2) < (1/2)1^2|_ Combining this 
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with the previous claim about the distribution of ID2I shows that the (now unconditional) 
probability of the event F that there is some x & D2 for which |\4| > is at least 
1 - (4 + /i2(e))-" - (l/2)^2n2 ^ _ ^ /i3(e))"" for < /12 and sufficiently large n. 

Finally, condition on F. Let x be an element of D2 with \Vx\ > v? jl and let D3 = T4. 
It suffices to show that with high probability there is some y E Ds such that y is infected 
at time t4, since we already know that y is dual infected at time t^. Recall that the only 

way that infection spreads between times and t4 is between neighboring vertices that 
intersect {rii + 1, . . . ,n} equally. This effectively breaks i3„ into 2"'~"'^ fibers on which 
propagation of the infection is independent and behaves like a Richardson's model of 
dimension rii. Each element of is in a different fiber, and uniformly for y E the 
probability that infection passes from y to the opposite corner y of the fiber in the time 
t4 — = 1 + e/7 tends to 1 for large n by Theorem 5.2 and in particular is eventually 
greater that 1/2. Thus the conditional probability given F that some y is infected at 
time t4 is at least 1 — (1/2)"^/^ for large enough n. 

Combining all the conditional probabilities and using Lemma 6.1 yields a probability 
of at least l-(4+/i3(e))-"-(l/2)"'/2 that i is infected by time ^3+^4 = 2 ln(4+2v^)+3+e. 
Pick h{e) < /13(e). Then when n is large enough, our bound is at least 1 — (4 + h{e))~'^ 
and the proof is finsihcd. □ 



6.2 The lower bound 

Theorem 6.4 For any e > 0, P(^(| ln(2 + v^) + In 2 - e) = B„) ^ as n ^ 00. 

Proof: First stochastically dominate Richardson's model at time t = ^ ln(2 + v^) — e/2 
by the corresponding value of a branching translation process, as in Section 5. Let x 
be any vertex with d{0,x) > n/2. Then in the notation of the proof of Theorem 5.3, 
P(a; e A{t)) < e'^'pix, t) < e"*((l - e'^*) /2)"/2((l + e'^*) /2)'^/2 = e'^*((l - e"^*) /4)"/2 = 
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[(e2*-e-2*)/4]". Now since e^* < v^S + 2, it follows that e^^-e'^* < ^/5 + 2-{V5-2) =4, 
so P {x E A(t)) uniformly in such a; as n — > oo. Since at least half the vertices of Bn 
satisfy d{6,x) > n/2, it follows that P{\A{ty\ > 2"-2) ^ 1. 

Now condition on > 2"'"^. The process {A{s) : s > is stochastically 

dominated by the process {Z{u) : u > t} for which Z{t) = A{t) and transitions from 
5" to 5" U {y} occur at rate n for all 5" and y ^ S. Now for each x G A{ty, P{x e 
Z{t + s)) = 1 — e""*, and furthermore these events are independent as x varies. Thus 
P{Z{t + s) = Bn\\A{tY\ > 2"-2) < (i_e-"^)2""'. Plugging in s = ln2-e/2 < ln(2-e/2) 
(for e < 2) and using 1 — e < e^*" gives 

PiZ{t + s)^Bn\ \A{ty\ > 2"-2) < [e-(2-/2)-"]2"-^ ^ e-(^)"/^ 
which tends to as n — > oo. The theorem now follows readily. □ 
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